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Abstract 

We study algebro-geometric (finite-gap) and eUiptic solutions of fully discretized KP or 2D Toda 
equations. In bilinear form they are Hirota's difi'erence equation for r-functions. Starting from a given 
algebraic curve, we express the r-function and the Baker-Akhiezer function in terms of the Riemann 
theta function. We show that the elliptic solutions, when the r-function is an elliptic polynomial, form 
a subclass of the general algebro-geometric solutions. We construct the algebraic curves of the elliptic 
solutions. The evolution of zeros of the elliptic solutions is governed by the discrete time generalization 
of the Ruijsenaars-Schneider many body system. The zeros obey equations which have the form of 
nested Bethe-Ansatz equations, known from integrable quantum field theories. We discuss the Lax 
representation and the action-angle-type variables for the many body system. We also discuss elliptic 
solutions to discrete analogues of KdV, sine-Gordon and ID Toda equations and describe the loci of the 
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1 Introduction 



Among a vast class ol solutions to classical non-linear integrable equations elliptic solutions play a special 
role. First, these solutions occupy a distinguished place among all algebro-geometric (also called finite-gap) 
solutions, i.e. solutions constructed out of a given algebraic curve. The general formulas in terms of Riemann 
theta-functions become much more effective - in this case the Riemann theta-function splits into a product of 
Weierstrass cr-functions associated to an elliptic curve. Second, there exists a remarkable connection between 
the motion of poles (zeros) of the elliptic solutions and certain integrable many body systems. 

The pole dynamics of elliptic solutions to the Korteweg - de Vries (KdV) equation and the Calogero- 
Moser system of particles were linked together in the paper § (see also §). It has been shown in Pl,|4|l 
that this relation becomes an isomorphism if one considers elliptic solutions of the Kadomtsev-Petviashvili 
(KP) equation. More recently, these results were generalized to elliptic solutions of the matrix KP and the 
matrix 2D Toda lattice equations (see js) and d), respectively). The dynamics of their poles obeys the spin 
generalization of the Ruijsenaars-Schneider (RS) model Q]. 

Let us recall some elements of the elliptic solutions for the standard example of the KP equation 5uyy = 
(4Mt -|- 6uUx — Uxxx)x for a function u = u(x, y, t). An elliptic solution in the variable x is given by 



N 



u{x,y,t) = const + 2^ p{x - x^{y,t)), (1.1) 

i=l 

where p{x) is the Weierstrass p-function. The self-consistency of this Ansatz is a manifestation of integra- 
bility. It has been shown in that the dynamics of poles as functions of y obeys the Calogero-Moser 

many body system with the Hamiltonian 



This system in it s tu rn is known to be integrable. There is an involutive set of conserved quantities H^^^ - 



the Hamiltonian (1^) and the total momentum are H^'^'> and H'-^\ The equations of motion are 

N 

d^Xi^A ^ p'{xi-Xj). (1.3) 

The ^-dynamics is described by H'^^\ 

The reduction to the KdV equation restricts the particles to the locus Lat in the phase space: 



-'N 



(k,.tO pi^0,^p'{x,~xj)^0^ (1.4) 



(here pi = dyXi). In spite of interesting developments, an analysis of the locus structure is far from to be 
completed. 

In this paper we extend these results to the fully discretized version of the KP equation or 2D Toda 
lattice. Being fully discretized they become the same equation. In bilinear form they are known as Hirota's 
bilinear difference equation (HBDE) |^ (see |^ for a review). This is a bilinear equation for a function 
r(/, m, n) (called T-function) of three variables: 

At(/ -|- 1, m, n)T{l, m -|- 1, 71 -|- 1) -|- fJ.T{l, m + 1, n)T(l + l,m,n + 1) 
+ i^T{l,m,n + l)T{l + l,m+l,n) = 0, (1.5) 

where A, fi, v are complex parameters and the three variables are not necessarily integer. In what follows we 
call them discrete times stressing the difference with continuous KP-flows. Let us introduce a lattice spacing 
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77 for one of the variables, say, I and denote x = rjl. By elliptic solutions (in the variable x) to this equation 
we mean the following Ansatz for the r-function: 



N 

t{1, m, n) = r'"^"(x) ^Y[a{x- x™'"), (1.6) 



where c^x) is the Weierstrass cr-function. We refer to the r.h.s. of (1.6) as elliptic polynomials in x. For 
brevity, we call solutions of this type elliptic though the r-function itself is not double-periodic. However, 
suitable ratios of these r-functions, for instance, 

rn,n( \ m+l,n( , \ 

A™'"(a;) = J'^i , (1.7) 

are already elliptic functions. 

In this paper we derive equations of motion for poles of A"^'^{x) or zeros of t™'"(x) for discrete times 
m, n and thus obtain a fully discretized Calogero-Moser many body problem. This appears to be the discrete 
time version of the Ruijsenaars-Schneidcr (RS) model proposed in the seminal paper [l0| . Remarkably, the 
discrete equations of motion have the form of Bethe equations of the hierarchical (nested) Bethe Ansatz. 
The discrete time runs over "levels" of the nested Bethe Ansatz. 

We also consider stationary reductions of HBDE. In this case the initial configuration of poles (zeros) 
is not arbitrary but constrained to a stable locus as in the continuous case ( |l.4| ). For the most important 
examples we give equations defining the loci. 

A renewed interest in soliton difference equations, and especially in their elliptic solutions is caused by 
the revealing classical integrable structures present in intcgrablc models of quantum field theory. It turns 



out that Hirota's equation (1.5) is the universal fusion rule for a family of quantum transfer matrices. Their 
eigenvalues (as functions of spectral parameters) obey a set of functional equations [ pT| which can be recast 
into the bilinear Hirota equation [|l2j (see also |l3j and for less technical reviews). Furthermore, it turned 
out that most of the ingredients of the Bethe Ansatz and the quantum inverse scattering method are hidden 
in the elliptic solutions of the entirely classical discrete time soliton equations p2]. In particular, the discrete 



dynamics of poles of A™'"(a;) or zeros of (1.6) has the form of Bethe Ansatz equations, where the discrete 
time runs over "nested" levels. 

The theory of elliptic solutions has direct applications to the algebraic Bethe Ansatz and to Baxter's 
T-Q-relation, which we plan to discuss elsewhere. 

Here we attempt to develop a systematic approach to the elliptic solutions of the integrable difference 
equations. The basic concept of the approach is the Baker- Akhiezer functions on algebraic curves. We prove 



that all solutions to HBDE of the form (1.6) are of the algebro-geometric type and present them in terms of 
Riemann theta functions. 

The plan of the paper is as follows. 

In Sect. 2 we describe general algebro-geometric (finite-gap) solutions to HBDE. We start from the Baker- 
Akhiezer function constructed from a complex algebraic curve of genus g with marked points. This function 
satisfies an overcomplete set of linear difference equations. Their consistency is equivalent to Hirota's equa- 
tion. In this way, one obtains a {Ag + l)-parametric family of quasiperiodic solutions to HBDE in terms of 
the Riemann theta-functions. Solitonic degenerations of these solutions are discussed in Sect. 2.4. 

Sect. 3 is devoted to elliptic solutions. They are shown to be a particular subclass of the algebro-geometric 
family of solutions of Sect. 2. We derive equations of motion for zeros of the r-function (the Bethe Ansatz 
equations) and their Lax representation. We discuss variables of the action-angle type and write down 
equations for the stable loci for the most important reductions of Hirota's equation. 
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2 Algebro-geometric solutions to Hirota's equation 



In this section we construct algebro-geometric solutions of Hirota's equation out of a given algebraic curve. 
The general method of constructing such solutions of integrable equations is standard. As soon as the bilinear 
equation can be represented as a compatibility condition for an overdetermined system of linear problems, 
the first step is to pass to common solutions ^ to the linear problems. Given a linear multi-dimensional 
difference operator with quasiperiodic coefficients, one associates with it a spectral curve defined by the 
generalized dispersion relations for quasimomenta of Bloch eigenfunctions of the linear operator. The Bloch 
solutions 5" are parametrized by points of this curve. Solutions to the initial non-linear equation are encoded 
in the analytical properties of ^ as a function on the curve. Spectral curves of general linear operators with 
quasiperiodic coefficients are transcendental and, therefore, intractable. However, soliton theory mostly deals 
with the inverse problem: to characterize specific operators whose spectral curves are algebraic curves of 
finite genus. Such operators are called algebro-geometric or finite-gap. Their coefficients yield solutions to 
HBDE, that we are going to study. 

Finite-gap multi-dimensional linear difference operators were constructed in the paper fT^ by one of the 
authors. We present the corresponding construction in a form adequate for our purposes. 

2.1 The Baker- Akhiezer function 

As usual, we begin with the axiomatization of analytical properties of Bloch solutions 'J. The Baker- Akhiezer 
function is an abstract version of the Bloch function. Since we solve the inverse problem, the primary objects 
are ^'-functions rather than linear operators. 

Let r be a smooth algebraic curve of genus g. We fix the following data related to the curve: 

- A finite set of marked points (punctures) Pa £T, a — 0,1, ... , M; 

- Local parameters in neighbourhoods of P^: vUaiPa) ~ 0; 

- A set of cuts Caf3 between the points Pa,P/3 for some pairs (it is implied that different cuts do 
not have common points other than their cndpoints at the punctures); 

- A set 13 of g (distinct) points 71, . . . , 7g g F. 

Further, we introduce the following complex parameters laf3 (times or flows): 

- To each cut Cafj is associated a complex number (it is convenient to assume that l^a = —lap)- 
Consider a linear space J-{1\ D) of functions "^{l; P), P G F, such that: 

1. The function ^'(^; P) as a function of the variable P G F is meromorphic outside the cuts and has at 
most simple poles at the points 7s; 

2. The boundary values \I>='=^("'') of this function at opposite sides of the cut Ca/3 satisfy the relation 



4'+'("'3)(Z;P) *-'("'')(/; P)e- 



(2.1) 
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In a neighbourhood of the point P^ 



a 



it has the form 




(2.2) 
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Note that if lap are integers, then ^' is a meromorphic function having simple poles at 7^ and having zeros 
or poles of orders \La\ at the points Pa- 

Any function 5" obeying conditions 1 - 3 is called a Baker- Akhiezer function. In our approach, these 
functions are central objects of the theory. In some cases (especially in the matrix generalizations of the 
theory) the notion of the dual Baker- Akhiezer function is also important (see e.g. Q, ||^). We omit its 
definition because it can be easily restored using [D , |^ . 

Let us prepare some notation. Fix a canonical basis of cycles ai,bi on F and denote the canonically 
normalized holomorphic differentials by duji, i — 1,2, . . . , g. We have 

dujj = 6ij , <j) dujj — Bij , 

where B is the period matrix. Given a period matrix B, the g-dimensional Riemann theta-function is defined 

by 

Q{X) e{X\B) = exp \ni{Bn, n) + 27ri(n, X) . 

Here X = {Xi, . . . , Xg) is a g-component vector. 

For each pair of points Pa,Pi3 G F, let dfl^°'^'> be the unique differential of the third kind holomorphic on 
F but having simple poles at the points Pa, Pp with residues —1 and 1 and zero a-periods. 

To write an explicit form of Baker- Akhiezer functions, let us choose one of the marked points, say Pq, 
andbyi*(P) = (Ai(F),... ,A<,(P)), 

MP) = r 

J Pa 

denote the Abel map. The Baker-Akhiezer function is given by the following theorem. 

Theorem 2.1 // the points 71, . . . , 7^ are in a general position (i.e. D is a non-special divisor), then J- {I; D) 
is a one- dimensional space generated by the function "^{l; P) G J' {I', D) 

^ e(l(p) + x(/) + z|B)9(z» / 

Q{A{P) + Z\B)Q[X[l) + Z\B) ' Jq. I 

Here Qo G T is an arbitrary point in the vicinity of Pq, Qq ^ Pq, belonging to the integration path Pq — > P 
in the Abel map. Further, 

Z = ^K-j2Mli). X{l) = Y.U^^PKp. (2.4) 
where K is the vector of Riemann's constants and the components of the vectors U^^f^^ are 

^^'""^^^'^ii dn^'''^ = MPp)~MPc.)- (2.5) 

The proof of theorems of this kind as well as the explicit formula for ^ in terms of Riemann theta-functions 
are standard in finite-gap theory (see e.g. flBl). The last equality in (2.5) follows from Riemann's relations. 



dui. 



Remark 2.1 Although the explicit formula ( |2.3| ) requires a fixed basis of cycles, the Baker- Akhiezer function 
is modular invariant. 



Remark 2.2 Since abelian integrals in (2.3) have logarithmic singularities at the punctures, one can define 
a single valued branch of ^! only after cutting the curve along Cap- 

Remark 2.3 The choice of the initial point of the Abel map is in fact not essential. It can be chosen not 
to be one of the punctures. In particular, it may be Qo, which slightly simplifies the theorem. However, our 
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choice simplifies the hnear equations for befow. Moreover, below we assume that the integration paths in 
the Abel maps A{Pa) go along the cuts Cqq. 

Remark 2.4 The theorem implies that any function from !F{1; D) has the form r(/)^(/; P), where r{l) is an 
arbitrary function of I but does not depend on P. It is convenient to choose it such that at the point Pq the 



first regular term in (2.2) equals 1. 

CoefRcients ^s"'' of the asymptotical behaviour of the '5 ( |2.2| ) can be expressed through the r-function 

T{l) = T{{l^p}) = e{X{l) + Z). (2.6) 

In particular, 

where Xa0 are Z-independent constants. Here and thereafter we skip unshifted arguments. 

Remark 2.5 If the graph of cuts includes a closed cycle, then a shift of variables la/i lap + 1 does not 
change the r-function but multiplies the ^-function by a cycle dependent constant. For instance, if the cycle 
consists of three links Cap, Cp^, C^a, then 

Ti^aP + Ij Ipf + Ij Ija + 1) = 'T{lap, Ipf, If a), 

(2.8) 

■^{lap + l,lp-f + I, If a + 1;^) = COllsm'{lap,lpf,l.ya;P)- 



This follows from ( p.lD 

In the sequel, we do not need the above construction in its full generality. For our purposes it is enough 
to consider the case of four punctures Pq, ■ • ■ , ^3 and a general graph of cuts as is in the figure. Cuts connect 
each pair of points. Any three links (not forming a cycle) give rise to a bilinear equation of the Hirota type. 
They have different forms, but are in fact equivalent due to ( ^.^ ). For further convenience we specify 

'01 = h, I02 = h, lo3 = I12 = h- (2-9) 




The general case of more punctures yields higher Hirota equations (i.e. the discretized KP or 2D Toda lattice 
hierarchies) . 
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2.2 Difference equations for the Baker- Akhiezer function 



The Baker- Akhiezer function 5'(Z;P) obeys certain Unear difference equations with respect to the variables 
laf3- Coefficients of the equations are fixed by the analytical properties of P) as a function of P G F. We 
restrict ourselves to the case of four punctures and use the notation introduced at the end of the previous 
subsection. The general case of more punctures can be treated in a similar way. 

Theorem 2.2 Let ^'(Z;P) be the Baker- Akhiezer function normalized so that ^g*^^ = 1. Then it satisfies the 
following linear difference equations: 



+ 1, l/i; P) - Ip + 1; P) + A^p[l^,lp)-^{l^,lp- P) = 



with 



(2.10) 
(2.11) 



for any a, P — 1,2,3, a ^ (3. 



The proof is standard in finite gap theory. Denote the l.h.s. of eq. (2.10) by 5'. This function has the same 
analytical properties as 4*. At the same time the leading term at the point Pi is zero: ^q^-* = for any la, 1/3- 
From the uniqueness of the Baker-Akhiezer function it follows that = 0. Eqs. (2.18), (2.2C) are proved in 
the same way. 



Remark 2.6 The dual Baker-Akhiezer function obeys difference equations obtained from eqs. (2.10), 
(2.18) by conjugating the difference operators in the right hand sides. 

The coefficient functions of eq. ( p. 10 ) are given by the leading nonsingular term ^q"-* of the Baker-Akhiezer 
function at the punctures. They can be found from eq. (|2.3|) and are expressed through the r- function (p^): 



T{la,lp)T{la + 1,1(3 + 1) 
T{la,ll3 + 1)t(/q + 1, 1(3) 



(2.12) 



The constants Aq,/3 are expressed through the constant terms r^"^'' in expansion of the abelian integrals 



[5^13 - 5^ a) log + r^"^^ +0{w^). 



as follows: 



= - exp {r, 



(0/3) _ ^(0/3)^ ^ 



(2.13) 



(2.14) 



It can be shown that A^q = —^ap for ol ^ [3 and 



Aa/3 



exp 



(2.15) 



for any cyclic permutation of {a/37} = {123}. The integration path goes from P~^ to the neighbourhood of 
Po along the cut Cg-y (in the opposite direction), passes through the point Qo then goes along the cut 
C'oq- 

Eqs. (^T0|) , (^Tl|) can be viewed as linear problems for the discretized KP equation. Choosing another 
triplet of variables, say ^1,^3, ^3 and using eq. ( |2.8[ ), i.e. 

T{h + 1, h;k + 1) = T{h,l2 + 1, h; h), (2.16) 
T{k,h,k-,l3)^e(Y,^iPc.)lc.+ {A{P2)-A{Pi))h + z] , (2.17) 
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one may rewrite eqs. (2.10), (2.11) in a form suitable for discretization of the 2D Toda lattice. In this case 
eq. ( ^.10 ) for a, (3 = 1,3 remains the same. Another linear equation. 



with 



^ihMP) - ^3 + i;P) + Bih,h)^ih - 1, h; P) = 0, 



B{h,k 



^^(/i-lJa) 



= -A 



12- 



T{h;h + l)T{h;h) 



follows from ( 2.10 ) for a, P — 1,2 as a result of the change of variables. The third equation, 

*(/3 + 1, [3 + 1; P) - i(/3, l3)-^ih, Is + l;P) = ^(^3 + 1, ^^3; P) - B{h, /"3)*(?3, ^"3; P) , 

with 

. J §^\h + lj3 + l)_ . T{h + l,l3 + l;k + l)T{h,k;h + l) 
^[h,^) — — rrr" — = ^ — — —-^12- 



^(1) (^3, [3 + 1) ■ '^^ riluh + 1; ^3 + iMh + 1, ^3; ^3 + 1) ' 

^^'^(;3 + l,r3) _ A12A32 t(Zi + 1,^3 + 1^3 + 1)^(^1, ^3^3) 



Al3 T{h,l3 + hhMh + 1, ^3; ^3 + 1) ' 



(2.18) 
(2.19) 

(2.20) 

(2.21) 
(2.22) 



is a linear combination of two other equations ( 2.1C| ) written in terms of the new variables. The constant 
prefactors in (2.21), ( 2.22| ) are derived using the reciprocity law for differentials of the third kind p^ . 



Alternatively, eqs. (2.18D, (2.20) can be proved in the same way as in Theorem 2.2. 



2.3 Bilinear equations for the r-function 

We have shown that the Baker- Akhiezer function satisfies an overdetermined system of linear equations. For 
compatibility of this system the coefficient functions must obey certain non-linear relations. In terms of the 
T-function all these relations have a bilinear form. 

Theorem 2.3 The r-function obeys the Hirota bilinear difference equation 

\2zT{h + 1, ^2, h)r{h,l2 + 1, ^3 + 1) - \i?.t{IiM + 1, h)T{h + 1, h, h + 1) 
+ Ai2r(?i,Z2,^3 + l)r(?i + l,;2 + l,/3)=0 (2.23) 



with the constants Xap defined in (2.14) 



Proof. First we examine the compatibility of eqs. (2.10) for (a/3) — (12) and (a/?) = (13). The variable h is 
common for both linear problems. Eqs. ( |2.10| ) can be rewritten as 



*(/„ + l;F) = Af^"^*(/„;F), a = 2, 3, 
where m}"^ is the difference operator in Zi: 

,,r(a) _ _ X T{h,lo)T{ll + 1,/a + 1) 

^ " ^^"r(Zi,Z„ + l)T(Zi + l,;„)- 



(2.24) 



(2.25) 



Eq. (2.24) has a family of linearly independent solutions parametrized by points of the curve F. Whence the 
compatibility is equivalent to commutativity of the operators 



[Mf\ Mf^] = 0, 



(2.26) 
(2.27) 
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which is a discrete zero curvature condition. Commuting the operators ( ^.25 ), we find after some algebra 
that this condition is equivalent to the relation 

Xl3T{ll,l2 + IJsMh + 1, h, 1-3 + I)- \i2t{IiM, h + + IM + 1, ^3) 

+ Hi{h;l2M)T{h + lMM)T{liM + lM + l)=Q. (2.28) 
where Hi is an arbitrary function such that i?i(/i + 1; ^3) = Hi{li; I2, Is). 

All this can be repeated for other two pairs of linear problems, i.e., for (a/?) = (21), (23) and (a/?) = 
(31), (32) in (2.10). This leads to bilinear relations similar to ( [2.28 ) for the same r-function but with different 
functions H2 and H3. To be consistent, all three bilinear relations must be identical. This determines 
Hi = — A23, etc, which proves the theorem. 



Remark 2.7 In the class of the algebro-geometric solutions the Hirota equation ( 2.23 ) is equivalent to Fay's 
trisecant identity [TtI] . 



The coefficients in eq. ( 2.2S ) may be hidden by the transformation 

^1^2 / \ \ hh 



Tih.hJs) 



Al3 

A12 



Al3 

A23 



T{h,l2,h) 



(2.29) 



bringing the Hirota equation into its canonical form: 

t{Ii + 1, h, h)r{llj2 + 1, ^3 + 1) - T{li,l2 + 1, h)T{ll + 1, /2, ^3 + 1) 

+ T{liMM + l)r{li + lM + l,h)^Q- 

The formula 

rP2 



T{liMM)=e^P\lil2 / dr!(°i'+?2?3 



P3 



(2.30) 



(2.31) 



thus provides the family of algebro-geometric solutions to eq. ( 2.3C ) . 

This family has 4g + 1 continuous parameters. Indeed, for 5 > 1 the solution depends on — 3 moduli of 
the curve, on g points "fi and on the 4 marked points Pa. The dependence on the choice of local parameters 
is not essential. 

In terms of the variables ^1,^3, ^3 the bilinear equations has the form 

Ai3r(/i,;3j3 + iMhJa + l-Js) - X23Tili,h;l3)T{li,h + l;h + 1) 

= Ai2r(;i + i,;3;r3 + i)T(Zi-i,;3 + i;/"3). (2.32) 

(Alternatively, this equation is a result of the compatibility of eqs. (2.1C), (2.18) an d (|2.20| ).) In contrast to 
eq. (2.23), the variable I2 is not shifted and skiped. The Hirota equation in the form ( 2.23| ) can be considered 
as a discrete KP, whereas the form ( 2.32| ) is a fully discretized 2D Toda lattice. Let us stress again that in 
the fully discretized setup the KP equation and the 2D Toda lattice become equivalent. 



2.4 Degenerate cases 

Degenerations of the curve F lead to important classes of s oluti ons. Among them are multi-soliton and 

rational solutions. Some examples of soliton solutions to eqs. ( 2.30 ), ( 2.32| ) were found by R. Hirota Here 



we outline the algebro-geometric construction of the general soliton solutions. 

Let us concentrate on multi-soliton solutions. In this case all N "handles" of the Riemann surface of 
genus TV become infinitely thin. In other words, the algebraic curve of genus N degenerates into the complex 
plane with a set of 2A^ marked points pi, g^: 

r {pi, Qt, « = 1, . . .,N}, 
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where pi, qi are the ends of the i-th handle. The Baker- Akhiezer function has the same value at each pair 
Pi, qi. The punctures Pa are replaced by points with local parameters Wa — z ~ Za- In this case the 
meromorphic differentials of the third kind have the form: 



[Zf3 - Za)dz 
{z - Zp){z - Za)' 

Let Fo{z) be a polynomial of degree N: 

N N 

Fo{z) = z^ + Y,b,z^-^=l[{z~j,). 

j=i i=i 

Its zeros ji will stand for poles of the Baker-Akhiezer function. 

Let us concentrate on the particular case when there are four punctures and three cuts from Zq to Za 
(a = 1,2,3) on the complex z-plane. Here Za are the z-coordinates of the marked points Pa'. Za = z{Pa). 
In this case the general definition of the Baker-Akhiezer function suggests the Ansatz: 

M,il;z)^^^^„il;z) (2.33) 

Fo{z) 



with 



and the polynomial 



a=l 



*oa;2)- n (^3^) ' l^ili,l2,h) (2.34) 



Z ~ Zq 



N 

F{l;z) = Y,aj{l)z^-^ (2.35) 

with yet undetermined l-dependent coefficients. These coefficients are determined by iV conditions 

^{hp,)^^{l;q,), Z = l,...,iV. (2.36) 

They are equivalent to the system of N linear equations for + 1 unknown coefficients aj in eq. ( 2.35| ): 

N 

J2K,,aj=0, i^l,...,N, j = 0,l,...,N, 

3=0 

where 

N-j N-j 

= IV^*o(''P») - |rrT^o(/;?.)- (2.37) 
FaiPi) Fo[qi) 

Solving the system of linear equations we represent the Baker-Akhiezer function in the form 

*(';^) = c(o||^*o(^;^), (2.38) 

where 

A{l;z) = det{Ki-i) 

ij 

and K is the {N -I- 1) x [N + l)-matrix with entries 

^Oj=Z^"^ K^J^K,J, i^l,...,N. 
The normalization factor c{l) is fixed by the asymptotics 

^{l- z)^{z- zo)-'^-'^-'^(l + 0{z - zo)) 
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near the point zq. It gives 



We also point out the identity 



A 0;zo -pr 



A(Z„;Zc<) = A(Za + l;2:o), a = 1,2, 3. 



The degeneration of the r-function ( 2.17 ) is 

t(/) = A(Z;zo). 



(2.39) 



(2.40) 



It satisfies the bihnear equation (2.23) with 

Aq/3 — 



{zo ~ Zq)(zo - Zji) 



The continuous parameters of the solution (|2.40| ) are IN points Pi,qi, N points ji and 4 points Zq. However, 
the r-function is invariant (up to an irrelevant constant) under the simultaneous fractional-linear transfor- 
mation z —^ (az + b)/ (cz + d) , ad — be = 1, oi all these parameters, so we are left with 3N + 1 parameters. 

In the case of rational degeneration the points pi and qi merge and the condition ( 2.36| ) becomes 



(DiVl>)(pO = 0, 



,N, 



(2.41) 



where are some differential operators in z with constant coefficients. This changes the form of the matrix 
K but eqs. (2.38), ( 2.39| ) remain the same. 



3 Elliptic solutions 

General finite-gap solutions of Hirota's equation for an arbitrary algebraic curve F with punctures Pa are 
quasi-periodic functions of all variables la- Below we construct a special important class of solutions for 
which the quantity (1/7) is doubly periodic in one of the variables. The r-function in this case is an elliptic 
polynomial (L6). We call them elliptic solutions. We show that the elliptic solutions also imply a spectral 
algebraic curve and are therefore a subclass of the algebro-geometric solutions of Sect. 2. 

Among all algebro-geometric solutions described in the previous section the elliptic solutions in one of 
the variables or their linear combination are characterized as follows. Let 



(3.1) 



be a vector field in the space of variables la/s and let U — X^q/s ^apU^°'^\ Let us transport the r-function 
(2.6) along the vector field dx and denote it by 



t{x) = e(Ux + X{l) + Z 



(3.2) 



Consider a set of algebraic curves with punctures Pa and cuts such that the vector U has the property: 

• There exist two constants llJi,uj2, Im (aj2/'^i) 7^ 0, such that 2ujaU, a — 1,2, belong to the lattice of 
periods of the holomorphic differentials on F: 



t{x + 2uja) = e' 



t{x) 



where ra, Sa are constants. 



11 



The T-function is then an eUiptic polynomial in the variable x. Due to the commcnsurability of U and the 
lattice of periods, solutions to the equation 

Q{Ux + f ) = (3.3) 

are Xi{Z) + 2JiUJi + 2J2(-U2, where Xi{Z) belong to the fundamental domain of the lattice generated by 2wi, 
2lu2 and Ji, J2 run over all integers. Therefore, 

N 

e(c7a; + Z) = e°i^+"^^'[](T(a;-a:,(f)) (3.4) 



with x-independent oi, 02 and N > g. 

The requirement of the ellipticity imposes 2g constraints on 4(7 + 1 parameters. Therefore the dimension 
of the family of elliptic solutions is 2g + 1 . 

Below we concentrate on a specific case where the ellipticity is imposed along Iqi = li by setting U ~ 
7^-1 [/(oi) and X = rjli where 77 is a complex constant. 

The logic of this section is opposite to the one of Sect. 2. Here we solve the direct problem and show that 
all solutions which are elliptic in any "direction" lap of the form (1.6) with a time- independent degree N ^ 
are of the finite gap type (3.2). 



3.1 Equations of motion for zeros of the r-function 

Let us show how to obtain the equations of motion for zeros of elliptic solutions by elementary methods. 
First of all, let us rename variables to emphasize the "direction" of ellipticity: 

li = rjx, I2 = n, = m, l^ = fh, (3.5) 

so that the r-function has the form ( |l.6| ). Let us now consider one of eqs. ( 2.10|) , say for a,(3 = 1, 2. In this 
equation all variables except x and m are parameters and we skip them wherever it does not cause confusion: 
a^r'" 2;™' r'"'"(a;) T™(a;) 

7n( \ m+l( I ^1 

V''"(x + r;) - Ai3 IZ, f-^^'^ix) = V^™+i(x) . (3.6) 

^ ' r'"+i(x)T"(a; + 77) v ^ w \ j 

Let us look for solutions of the form 

^ ^ ' r™(a;) 
with some function p™(x). Eq. ( |3.6| ) then reads 

T 



™+i(x)p™(a; + 77) - Ai3r'"+i(:E + ri)p^(x) = t™(x -f- ^)p'^^\x). (3.7) 



We are interested in the case when r™(a;) is an elliptic polynomial in x for any m: 



r'"(x) = J|a(x-a;f). (3.8) 



The "equations of motion" for its roots a;™ in the "discrete time" m can be easily obtained from eq. (3.7) in 

1 

Ai3T"+^(a;f +^ + ?7)p"(:rf +^) ^ t"(x7+^ + r7)p'"+^(a;™+^). 



the following way. Substituting x = x™" , x ~ x^ — rj, x — xj^ ~ 77, we get the relations 



T 



^(x7+i - r7)p'"(x7+^) = r"(a;™+i)p"+i(:r7+i - 77), 



r™+i« - 77)P"«) = Xi3T"-^+\xf)p"\xf - ?7), (3.9) 
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respectively. Combining these relations, we eliminate p and obtain a system of equations for the roots x'^ 



a{x^ 



N 



')aixf ~x^ + r,)aix"^ -xl^+' -rj) 



-ri)a{x'; 



'n)a{x" 



-1. 



(3.10) 



Note that these equations involve only one discrete variable m while the direct substitution of the Ansatz 



(3.8) into the non-linear equation (2.30) would give a system of equations involving two discrete variables. 
It is not easy to see that they in fact decouple. The decoupling becomes transparent if one starts with the 
auxiliary linear problem (3.6). 



Similar equations hold for the discrete I3 
equations for the fh dependence of Xi : 



m-dynamics. From the linear problem (2.18) we obtain the 



N ( fi 



^)a{x^ 



m+l\ 



„m— 1 



?7)cr(a;^ 



v) 



-1. 



(3.11) 



To avoid confusion, let us stress that the Xi depend on all discrete times. However, the variables are separated 
in the equations of motions. 

Remark 3.1 Another choice of the "direction" of ellipticity gives rise to different equations of motion. 
To illustrate this, let us require t{1) to be an elliptic polynomial in the direction orthogonal to s = I3 + I3 
and a = h + I3 — I3, so that the zeros of the r-function depend on a and s. Then the vector field in (3.1) is 
dx = — ^c?;., + ^9^3. In terms of T°-''{x) = t{Ii,13,13) bihnear equation ( 2.32| ) reads: 



Ai2r'^'^(a; + ?7)r"'^(x - ry) + X23T''''+Hx)T''''~\x) = 
and the zeros x^'' of r'^'*(a;) obey the system of coupled equations: 



'{x)T''-^''{x) 



v) 



-7;)r°.''+i(a:°'' -77) 

■7j)T-'^+\xr+V) 



ya-i,.+i(^°:«)ya+i,.(.^^.« + 7?)T'^.-i(a;- - 77) 

These equations are more complicated than ( 3.10| ). In contrast to the previous case, evolutions in a and s 
are not separated. 



3.2 Double-Bloch solutions to the linear problems 



In order to further examine the elliptic solutions, we need the notion of double-Bloch functions. A mcro- 
morphic function f{x) is said to be double-Bloch if it enjoys the following monodromy properties: 



f{x + 2uJa)=Baf{x), 



a = 1,2. 



(3.12) 



The complex numbers Ba are called Block multipliers. A non-trivial double-Bloch function can be represented 
as a linear combination of elementary ones: 



N 



/(x)=^c,$(a;-a;„C)fc"/^ 



(3.13) 



where 



a{( + X + 7]) 
a{( + T])cr{x) 



^(C - V) 



x/(2vi) 



and complex parameters C and k are related to the Bloch multipliers by the formulas 

V(C-ry)X-/" 



= fc2-./"exp(2C(u;„)(C + 77)) 



a{C + ri) 



(3.14) 



(3.15) 
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— a' {x) / (t{x) is the Weierstrass C- function). 

Let us point out some properties of the function $(a;,C). Considered as a function of C, ^{xX) is 
double-periodic: 

<i>(x,C + 2c^a)=<i>(x,C). 

For general values of x one can define a single- valued branch of <i>(x, C) by cutting the elliptic curve between 
the points C = ±77. In the fundamental domain of the lattice generated by 2LJa the function Q) has a 
unique pole at the point a; = 0: 

$(x,c) = - + o(i). 

X 

In the next subsection we need the identity: 

<^{x, z)$(y, z) = <^{x + y, z){C,{x) + Civ) C,{z + ii) - C,{x + y + z + ry)) (3.16) 

which is equivalent to the well known 3-term bilinear functional equation for the u-function. 

Recall the notion of equivalent Bloch multipliers The "gauge transformation" f{x) f{x) = f{x)e^^ 
{b is an arbitrary constant) does not change the poles of any function and transforms a double-Bloch function 
into another double-Bloch function. If Ba are Bloch multipliers for /, then the Bloch multipliers for / are 
Bi = Bie^^'^^, B2 = -626^'"^^. Two such pairs of Bloch multipliers Ba and Ba are said to be equivalent. (In 
other words, they are equivalent if the product B'^^B2'^^ is the same for both pairs.) 

This definition implies that any double-Bloch function can be represented as a ratio of two elliptic 
polynomials of the same degree multiplied by an exponential function and a constant: 

i—l ^ ^ 



The Bloch multipliers are 



N 



Ba = (fc')'-"/" exp I 2C(c.,) Y^{x, - y,) 

3 = i 



Eqs. (3.13) represents a Bloch function by its poles and residues, whereas eq. (3.17) represents a Bloch 
function by its poles and zeros. 



3.3 The Lax representation 



The coefficients in eq. ( |3.6D are elliptic functions, i.e. double-periodic with periods 2uJa- Therefore, the 
equation has double-Bloch solutions . Similarly to the case of the Calogero-Moser model and its spin 
generalizations the dynamics of poles of the elliptic coefficient in the linear problem is determined by the 
fact that equation (3.6) has an infinite number of double-Bloch solutions. 

In what follows we always assume that the poles are in a generic position, i.e. — a;™ 7^ 0, ±77 and 
X™ — x™^^ 0, ±77 for any pair i ^ j. Exceptional cases are also of interest but must be treated separately. 



Theorem 3.1 LetT™{x) be an elliptic polynomial of degree N . Equation (3_^) has N linearly independent 
double-Bloch solutions with simple poles at the point s x™ and equivalent Bloch multipliers if and only if zeros 
X™ of the T-function satisfy "equations of motion" (SAL). 



Theorem 3.2 // eq. (S.t ) has N linearly independent double-Bloch solutions with equivalent Bloch multipli- 
ers, then it has infinite number of them. All these solutions have the form 



N 



^^{x)=Y^c,{mX,m{x-xTX)k 



x/r] 



(3.18) 



i=l 



('$(x, C) is defined in (3.14))- The set of corresponding pairs (C, z) is parametrized by points of an algebraic 
curve. 
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These theorems are proved by the same argmiients as in [gj. Here we present the main steps. 

N Unearly independent double-Bloch solutions with equivalent Bloch muhiphers may be written in the 
form with some values of the parameters Cr, kr, s — 1, . . . , N. Equivalence of the multipliers implies 

that the can be chosen to be equal C,r = C- 



Let us substitute the function ^"^{x) of the form (3.18) with this particular value of C into eq. (p.q). 
Since any double-Bloch function (except equivalent to a constant) has at least one pole, it follows that the 
equation is satisfied if its left hand side has zero residues at the points x = x™' — rj and x = x^~^^ . The 
cancelation of poles at these points gives the conditions 



N 



kc,{m, C, k) - f,{m) C, - - 77, C) = , 



N 



c,{m + 1, C, k) = g,{m) ^ c,(m, C, k^xT^' - x^, () , 



where 



fi{m) = A 
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7])ct(x^ 



X 







y,rn+l 



- ri)a{x. 



m+l 
s 

,m+l 



v) 



x: 



m+l 



m+l 



v) 



Introducing a vector C{m) with components Ci{m, C, z) we can rewrite these conditions in the form 

{C{m)-kI)C{m)=Q, (3.23) 

C{m + l) = M{m)C{m), (3.24) 
where / is the unit matrix. The matrix elements of C{m) and Ai{m) are: 

C,,{m) = MmMxT - ~ 77,0, (3.25) 

A^,,(777) = .g,(777)$(x™+i - ,C). (3.26) 
The compatibility condition of (3. IE) and ( 3.20| ), 

li{m^\)M(m)= M(m)L(ra) (3.27) 

3T^) 



has a form of the discrete Lax equation. The Lax equation ( 3.27 ) appeared in ref. |10| , where eqs. 
were proposed as a time discretezation of the RS model. 



3.19) 



3.20) 



3.21) 



3.22) 



Lemma 3.1 For matrices C and M defined by {3.2L ), (3.2t ) the discrete Lax equation (3.21) is equivalent 
to the equations (3.1C). 

The proof is along the lines of ref. ||lfl . We have 



(x, C)iMim)C{m) - C{m + l)X(77i)), 



/j(777+ l)^.g,(777)$(5 



.m+l 



m+l 



7/,C)$(x 



m+l 



+ 9^{m) J2 MmMxT+' ~ x:\ C)$(a;r - xf - 77, C). 

The coefficient in front of the leading singularity at = —rj is proportional to 

9s{m) + gi{m) ^ /s(r77). 



(3.28) 
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On the other hand, 



^(/,(m)-5,(m)) =0 



(because this is the sum of residues of the ehiptic coefficient in eq. (3.6) ). Therefore, 

/i(m + 1) = 5j(to), i = l,...,N. 
These equations are equivalent to ( 3.10| ). 



To show that ( |3.28 ) is identically zero provided eqs. ( 3.30 ) hold, we use the identity ( 3.16 ) 



(3.29) 



(3.30) 



v)) 



v) + Cix. 



m+l 



'))■ 



(3.31) 



Noting that G is proportional to the sum of residues of the elliptic function 



N 



[axT+'-v-x)+ax-xf)]'[[ 



a{x - x'^')a{x - x'l'+' + r/) 
^Jl a{x ~x'^ + ri)a{x - a;f +i) 



at the points x = and x = — 77, we conclude that G = 0. 



It was already proved that eq. (3.6) has N linearly independent solutions if eqs. (3.10) or the Lax equation 



(3.27) hold for some value of the spectral parameter C,. It then follows from Lemma 3.1 that the Lax equation 
holds for any value of C,. Therefore, for each C, there exists a double-Bloch solution given by ( 3.18| ), where 
the Ci are components of the common solution to (3.23), (3.24). 



Theorem 3.3 All elliptic solutions of eq. (2.2c ) of the form (3.1 ) are of the algebro-geometric type and x^ 
are given implicitly by the equation 



e 



(^7j-^A{P,)xT+mA{P3) 



0, 



(3.32) 



where the Riemann theta-function corresponds to the algebraic curve F defined by the characteristic equation 

N 



R{k, C) = det(/:(TO) + kl) ^k^ + J2 niOk'^^' = . 

i=l 

The matrix C is given by eqs. (3.21 ), ( 3.2f\ ) and the coefficients ri(Q have the form 



(3.33) 



_ a(C + ^)^'-^)/V(C-(^-l)^) , 



where li are integrals of motion. The characteristic equation (3. 35) and Ij are functions of x"^ and a;™^^ 
but stay the same for all m. The spectral curve F dete rmin ed by eq. ( 3. 3^ ) is an algebraic curve realized as 
a ramified covering of the elliptic curve. The function ( 3.11 ) is the Baker- Akhiezer function on F. 



We call the spectral curve F defined in Theorem 3.3 the Ruijsenaars-Schneider (RS) curve. The RS curve is 
identical to the spectral curve for the continuous time RS model studied in rcf. Q . The proof of the theorem 
is omitted. It is as in ref . H . 



The matrix C is defined by fixing a;™° and x 



mo + l 



These Cauchy data uniquely define the RS curve F, 
the vectors A{Pa), a = 1,2 and Z in eq. (3.32). The curve and the vectors A{Pa), a = 1, 2 do not depend 
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on the choice of mg. They are action-type variables. The vector Z depends Hnearly on this choice and its 
components are thus angle- type variables. 

Remark 3.2 The discrete time dynamics defined in Theorem 3.3 is time-reversible, i.e. the Cauchy data 
^mo ^ 2:™o+i completely determine the dynamics for both time directions up to permutation of the "particles" . 
The £-A^-pair for the backward time motion is obtained from the difference equations for the dual Baker- 
Akhiezer function (see Remark 2.6) with an Ansatz similar to ( 3.1§| ). An alternative way to derive equations 
of motion ( p. 10 ) is to require the spectral equation (3.23) to be identical to the similar equation for the 
backward time motion. 



Remark 3.3 The form of equations for the dynamics in ^2 = n is identical to the equations ( |3.10 ) of the 
dynamics in to = ^3. The Cauchy data for TO-dynamics, i.e., values of at m = and to = 1 completely 
determine an evolution and Cauchy data in n (as well as all other flows). Comparing /^-operators for each 
flow, one finds relations between the Cauchy data: 



N I 0.0 



,.1,0 



0,0 



^0.1 



ry) 



J. X / 0,0 u,i\ / u,u i,u \ 
s=l -^s' -^V) 

Similar connections exist for the initial data of the m-flow. 



1,0 



A12 
Al3 ' 



z = l,2. 



,N. 



(3.34) 



3.4 Loci equations 

The values of x^, x} may be arbitrary if no other reduction (apart from the elliptic one) is imposed. If 
there is an additional reduction, then the x] are constrained to belong to a submanifold of C^^, the 



reduction locus. An example of such a locus in the continuous setup is the KP— i>KdV locus (1.4). Here we 
present equations defining the loci for three important reductions of Hirota's difference equation. In these 
cases spectral curves of algebro-geometric solutions are hyperelliptic. As before, a;°, xj are assumed to be in 
generic position. 

A) Discrete KdV equation JT^/. Discrete KdV equation appears as the reduction 

t(/i, ^2 -I- 1, /a + 1) = T{li,l2, k), X = -qli 

of the general 3-dimensional Hirota equation ( ^.23 ). In the notation (3.5) the equation is 



A23r"(x + 77)T"(a;) - Ai3T'"-i(x)r"+i(x + 77) + Ai2T"+i(x)r"-i(a; + 77) . (3.35) 

For algebro-geometric solutions the reduction means that {7^"^-' -1- {/'"■^•' belongs to the lattice of periods. 
Therefore, the function 

'"''(df]("2)+df^(°3)) ] (3.36) 



e(F) = cxp 



is meromorphic on the curve F. From the definition of the abelian integrals it follows that this function has 
a double pole at Pq and simple zeros at P2 and P3. Outside these points e{P) is holomorphic and is not 
zero. The existence of such a function means that the spectral curve is hyperelliptic. A hyperelliptic curve 
of genus g can be defined by the equation 

2g+l 

n(^-^^)- (3.37) 

This is a two-fold covering of the complex plane of the variable e. The projection of F onto the e-plane 
defines £ as a meromorphic function on F. This function has a double pole on F at the branch point Poo 
(above e = 00) and two simple zeros at the points Pq^^ (above e = 0). The identification of this notation 
for the punctures with our previous ones is 

Poo = Po, Pt^ = P2, Po'+' = Pi- 



ll 



The branch points Si in (3.37) may not be arbitrary since the curve should simultaneously be of the RS 
type. Correspondingly, the Cauchy data x^, xj with respect to the Za-flow obey certain constraints. Using 
the equations of motion (3.10) and (3.34), we obtain a system of 27V coupled equations on allowed values of 
x^, x] ("equilibrium locus"): 



N 

n 



a{xl 



f])a{xl 



v) 



N 

n 



cO + T^)aix'. 



v) 



V) 



?7)cr(a;0 



V) 



Al3 

Ai2 
Al3 



(3.38) 
(3.39) 



for i = 1, 2, . . . , iV. With the help of eq. ( 2.15 ) the right hand sides can be expressed through abelian integrals. 
The relation between the number of zeros N and genus g of the spectral curve is a subtle question. We do 
not discuss it here. 

Each of the systems ( 3.3^ ), ( |3.39| ) has the form of Bethe equations for an iV-site spin chain of spin 1 at 
each site. One may treat (fo r instance) as arbitrary input parameters while Bethe's quasimomenta xl are 
to be determined by eqs. ( 3.38 ). However, the system of locus equations ( 3.38 ), ( 3.39 ) determines x^ and xj 
simultaneously. 



Xi , X- 



12" 

Xi -\- €Xi 4" ^6 X j , Xl 



In the continuous time limit we set xf 
A13/A12 —(1 + Ce) with a constant C, we get from (B.38), ( 3.39| ) 

N 

^ Xk [C,{xi - Xk -1]) - C{xi -Xk+ 77)] C, 



dtXi 



fc=i 

N 



ik [p{xi -Xk-rf) - p(a 



+ ri)] = 0. 



0. Assuming 
(3.40) 
(3.41) 



fc=i 



In the leading order in e the systems (|3.38 ), ( 3.39| ) coincide with each other and yield the first system ( 3.40| ) 
while the second one (3.41) follows from the higher order terms. These equations define the equilibrium 
locus for the Ruijsenaars-Schneider system of particles. 

B) ID Toda chain in discrete time The reduction condition in this case is: 

T{lij2,h + 1,^3 + 1) = T^h^hJaJs)- 
The discrete time ID Toda chain in the bilinear form reads 



A 



13^ 



(xy+'ix) - Xi2t"'-'{x + 77)t"+^(x - 77) = A23(r"(2:)) 



(3.42) 



where we have excluded I3 and have passed to the notation of Example A). In this case U^-^^^ + [7^^^^ belongs 
to the lattice of periods. The corresponding curve is given by eq. (3.37) with a polynomial of even degree in 
the r.h.s. 



The Cauchy data x^, x} with respect to the l^-How satisfy the locus equations: 

N 

n 

s=l 
N 

n 



a{x} ~xl+r])a^ixl 


-x'i) 


A12 


'^{xj -xl~ i])a'^{x} ~ 


x° + v) 


Al3 ' 


aix°^ - x° + 7^)a\x° - 


xl ~ v) 


A13 


cr(x" — x^ — r/)a'^{Xi 


-xl) 


A12 



(3.43) 
(3.44) 



The continuous time limit can be taken similar to the way of the previous example. In this case, however, 
we have to assume A13/A12 Ce^^ as e — > with a constant C. We get 



N 



n 



ajx^ - Xk+ i])(j{xi - Xfc - 7?) 

0-2 (xi - Xk) 



= Cxh 



(3.45) 
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N 



(xi + Xk) [C,{xi - Xk -1]) + C{xi -Xk+Tl) - 2C{xi - Xk)] = 0. (3.46) 



These equations follow also from eqs. (4.58), (4.59) of the paper They define the stable locus for the RS 



system with respect to another flow than in eqs. (3.40), (3.41). 
C) Discrete sine-Gordon equatior^. The reduction condition is 

T{li,h, h + 1, ^"3) = t(/i + 2, 12, k, + 1), 
so, passing to the same independent variables as in the previous examples, we get the equation 

Ai3r'"(:r - 77)T™(a; + v) ~ A23r'"-i(x + 7/)r"+i(x - ??) = X,2{t'"{x))^ . (3.47) 

Now it is the vector [7'^°^^ + ^Yied belongs to the lattice of periods. The continuous SG equation is 

reproduced in the limit Pq, P2 Pi- The spectral curves are again hyperelliptic. 

The locus equations for the Cauchy data x^, xj with respect to the m-flow are 

o-(a:° -xl)a{x° -xl -2?]) _ A12 , . 

aHxO-xi-,) - xTs' ^'-''^ 

^ a(x[-x>(^[-x2 + 2^ _ A12 

n a^xl^-xo + r^) -Ai3- ^''"'^ 

Note that the structure of these equations is different compared to the previous examples. 

3.5 Remarks on trigonometric solutions 

Trigonometric solutions are degenerations of the elliptic solutions when one of the periods tends to infinity. 
They form a particular subfamily in the variety of soliton solutions described in Sect. 2.4. The trigonometric 
solutions admit a very explicit description in terms of the data defining the singular curve. 

Let us set the period to be 27r: 

T'"(a; + 27r) :=T™(a;), 



so an elliptic polynomial becomes a Laurent polynomial in e'^. The Bethe-like equations on motion (3.1C) 
preserve its form, but the Weierstrass function (t{x) is replaced by sin a;. 



It follows from the periodicity that the function "^q (2.34) obeys 

j = l,2,...,iV, 



^o(?7X + 27r?7,m;pj-) _ -^gi^qx^m-^Pj] 



■^oivx + 2771], m;qj) ^'0(772:, to; q^) 
or, explicitly 

P^^q^^-.,j,^ (3.50) 
Pj - Zo Qj - Zo 

where Jj are integer numbers. This condition restricts admissible singular curves (Riemann spheres with N 
double points). Minimal Laurent polynomials correspond to the choice Jj = ±1. This gives TV conditions 
for Pj, Qj, so the number of continuous parameters in the trigonometric case is 2N + 1 - the same as for the 
non-degenerate curves of genus N . 

Here we do not discuss the trigonometric degeneration of the >C-A4-pair and dependence of of pj and qj 
on initial data for Ruijsenaars-Schneider particles. This can be done along the lines of the paper 

^This version of the discrete SG equation is diflferent from the one considered in |^], [po| . The latter is closer to a special 
degeneracy of the discrete KdV at A12 — » A13. 
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The trigonometric loci can be characterized alternatively by imposing a relation on pj and qj in addition 



to (3.50) 



E{p,) = E{q,), p,^q,, i = l,2,...,N. (3.51) 
For the examples A)-C) of Sect. 3.4 the functions E{z) are 

A) E{z) - - - 

B) ij;(z)^ (3-52) 

C) E{z) 





(z - zo 


? 




- ^2)(^ 


-Z3) 










-z^){z 


-Z2) 






- zq) 



Conditions ( 3.5C| ), ( [3.51 ) leave us with a discrete set of admissible pairs Pi,qi. The continuous parameters 



7s give then an implicit paramctrization of the loci. 

4 Conclusion 

We have shown that the main body of finite-gap theory and the theory of elliptic solutions to nonlinear 
integrable equations is also applicable to finite-difference (discrete) integrable equations. Discrete equations 
includes the continuous theory as the result of a limiting procedure. Furthermore, discrete equations reveal 
some symmetries lost in the continuum limit. We have shown that all elliptic solutions with a constant 
number of zeros in the evolution (compare to [p^), are of the algebro-geometric type. Moreover, their 
algebraic curves are spectral curves for L-operators of the Ruijsenaars-Schneider model. Each point of this 
curve gives rise to discrete time dynamics of zeros of the r-function. 

The structure of equilibrium loci equations of reductions of Hirota's equation (analogues of the known 
KdV- locus ( |1.4| )) is expected to be richer than in the continuous case and requires further study. It would 
be very interesting to extend the algebro-geometric approach to elliptic solitons of KdV of ref. to the 
difference case as well as to understand difference elliptic solitons in terms of the Weierstrass reduction theory 

H- 

To the two main motivations pointed out at the beginning of the paper we can now add yet another one: 
an intriguing intimate connection between the elliptic solutions to soliton equations and quantum integrable 
models solved by the Bethe Ansatz. In our opinion, the very fact that the zeros dynamics and equilibrium 
loci are described by Bethe-like equations is remarkable and suggests hidden parallels between quantum and 
classical integrable equations. 



Acknowledgments 

We thanks Ovidiu Lipan for discussions and for his interest in the subject and J.Talstra for help. The work 
of I.K. was supported by RFBR grant 95-01-00751. The work of A.Z. was supported in part by RFBR 
grant 97-02-19085, by ISTC grant 015 and also by NSF grant DMR-9509533. A.Z. is grateful to Professor 
R. Seller for the kind hospitality at the Technische Universitat Berlin where a part of this work was done. 
P.W was supported by NSF grant DMR-9509533. P.W and A.Z thank Institute for Theoretical Physics at 
Santa Barbara for its hospitality in April 1997 where this paper has been completed. 



References 

[1] H.Airault, H.McKean and J.Moser, Rational and elliptic solutions of the KdV equation and related 
many-body problem, Comm. Pure and Appl. Math. 30 (1977), 95-125. 



20 



D.V.Chudnovsky and G.V.Chudnovsky, Pole expansions of non-linear partial differential equations, 
Nuovo Cimento 40B (1977), 339-350. 

I.M.Krichever, On rational solutions of Kadomtsev-Petviashvilii equation and integrable systems of N 
particles on line, Funct. Anal i Pril. 12:1 (1978), 76-78. 

I.M.Krichever, Elliptic solutions of Kadomtsev-Petviashvilii equation and integrable systems of particles, 
Func. Anal. Appl. 14:4 (1980), 282-290. 

I.Krichever, O.Babelon, E.Billey and M. Talon, Spin generalization of the Calogero-Moser system and 
the Matrix KP equation, Preprint LPTHE 94/42. 

I.M.Krichever and A.V.Zabrodin, Spin generalization of the Ruijsenaars-Schneider model, non-ab elian 



2D Toda ch ain and representations of Sklyanin algebra, Uspekhi Mat. Nauk, 50:6 (1995) 3-56, hep- 



th/9505039 . 



S.N.M.Ruijsenaars and H. Schneider, A new class of integrable systems and its relation to solitons, Ann. 
Phys. (NY) 170 (1986), 370-405. 

R.Hirota, Nonlinear partial difference equations II; Discrete time Toda equations, Journ. Phys. Soc. 
Japan 43 (1977) 2074-2078; Discrete analogue of a generalized Toda equation, Journ. Phys. Soc. Japan 
50 (1981) 3785-3791; 



A.Zabrodin, A survey of Hirota's difference equations, preprint ITEP-TII-10/97, ^olv-int/9704001 



F.Nijhof, O.Ragnisco and V.Kuznetsov, Integrable time-discretization of the Ruijsenaars-Schneider 
model, Commun. Math. Phys. 176 (1996) 681-700. 

P.P.Kulish and N.Yu.Reshetikhin, On GL^-invariant solutions of the Yang-Baxter equation and associ- 
ated quantum systems. Zap. Nauchn. Sem. LOMI 120 (1982) 92-121, Engl, transl.: J. Soviet Math. 34 
(1986) 1948-1971; A.Klumper and P.Pearce, Conformal weights of RSOS lattice models and their fusion 
hierarchies, Physica A183 (1992) 304-350; A.Kuniba, T.Nakanishi and J.Suzuki, Functional relations 
in solvable lattice models, I: Functional relations and representation theory, II: Applications, Int. Journ. 
Mod. Phys. A9 (1994) 5215-5312. 

I.Krichever, O.Lipan, P.Wiegmann and A. Zabrodin, Quan tum integrable models and discrete classical 
Hirota equations, preprint ESI 330 (1996), hep-th/960408C , to appear in Commun. Math. Phys. 

A. Zabrodin, Dis crete Hirota's equation in quantum integrable models, preprint ITEP-TH-44/96 (1996), 
tiep-th/9610039| . 

P.Wiegmann, Bethe Ansatz and Classical Hirota Equation Int. J. Mod. Phys. B, 11 (1997) 75-89. 

I.M. Krichever, Algebraic curves and non-linear difference equation, Uspekhi Mat. Nauk 33 (1978), n 
4, 215-216. 

B. Dubrovin, V. Matveev and S. Novikov, Non-linear equations of Korteweg-de Vries type, finite zone 
linear operators and Abelian varieties, Uspekhi Mat. Nauk 31:1 (1976) 55-136; I.M. Krichever, Nonlinear 
equations and elliptic curves. Modern problems in mathematics, Itogi nauki i tekhniki, VINITI AN USSR 
23 (1983). 

J. Fay, Theta functions on Riemann surfaces. Springer Lecture Notes in Math., 352 (1973). 

R.Hirota, Nonlinear partial difference equations I, Journ. Phys. Soc. Japan 43 (1977) 1424-1433. 

R.Hirota, Nonlinear partial difference equations III; Discrete sine-Gordon equation, Journ. Phys. Soc. 
Japan 43 (1977) 2079-2086. 

[20] L.D.Faddeev and A.Yu.Volkov, Quantum inverse scattering method on a spacetime lattice, Teor. Mat. 
Fiz. 92 (1992) 207-214 (in Russian); 



21 



[21] A.Trcibich and J.-L.Vordicr. Solitons elliptAques, Grothcndicck Festschrift, ed.: P.Carticr ct al, Progress 
in Math. 88, Birkhauser, Boston, 1990; A.Treibich, Tangential polynomials and elliptic solitons, Duke 
Math. Journ. 59 (1989) 611-627. 

[22] E.Belokolos and V.Enolskii, Verdier elliptic solitons and Weierstrass reduction theory, Funct. Anal, 
i Pril. 23:1 (1989) 57-58; E.Belokolos, A.Bobenko, V.Enolskii, A.Its and V.Matveev, Algebraic- 
geometrical approach to nonlinear integrable equations, Berlin: Springer, 1994. 



22 



